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^ Abstract. Utilizing asymptotic Pade-approximant methods, we estimate the three-loop order 
MS coefficients of a^[log(/i^/m^(/i)]'"' terms [k — {0, 1,2,3}] within the h ui^Vg decay rate. 

<^ • Except for the coefficient of the k = term, all other coefficients may also be obtained via 

. renormalization-group (RG) methods. The relative errors in asymptotic Pade-approximant esti- 

1^ ! mates of the k = {1, 2, 3} terms are found to be 5.1% or less, thereby providing an estimate of 

00 I the theoretical uncertainty in the asymptotic Pade-approximant estimate of the i?G-inaccessible 

' fc = term. By a judicious choice in the renormalization scale parameter, we are able to extract 

' \Vub\ from the inclusive decay rate to within ±7% theoretical uncertainty. 

-g^l INTRODUCTION 

The inclusive semileptonic B Xui^Vi decay rate is sensitive only to a single 
^ ■ mass parameter, the 6-quark mass rrib, since the w-quark, charged lepton, and 

>■ . associated antineutrino are all very light compared to nib. This property makes 

^ I the calculated decay rate particularly attractive for the extraction of the CKM 

^ I mixing parameter \Vub\- QCD corrections to the rate have recently been calculated 

to two-loop order within an MS scheme [1] . If the renormalization scale fi is chosen 
to be nibirnb) = 4.2 GeV^ [2], i.e., if /i = mb{fJ,), then the two loop rate is given by 

T{B XJ-Ve) = Kmlimb) [1 + 0.30 + 0.14] , K = Gl\VubWl92TT\ (1) 

Determinations of \Vub\ based upon comparison of experimental results to this ex- 
pression are necessarily subject to theoretical uncertainties associated with the 
premature truncation of the slowly convergent perturbative series in (1). In the 
present work, we demonstrate how asymptotic Pade-approximant methods in con- 
junction with renormalization-group (RG-) methods can dramatically reduce such 
theoretical uncertainties . . . 



1) by providing an estimate of three loop contributions to the decay rate that 
can be tested against RG-accessible three-loop coefficients in the rate, and 



2) by demonstrating that the physically-motivated choice of scale parameter /i 
which minimizes scale sensitivity [3] is quite close to that which minimizes 
the size of the four-loop term, and consequently, the truncation error for an 
asymptotic series. 



DETERMINATION OF RG- ACCESSIBLE 
THREE-LOOP COEFFICIENTS 



The series portion S[x,L\ of the perturbative MS expression for the inclusive 
semileptonic rate. 



r (5 ^ XJ-Ve) = K [m^V)]''5[a;(/.),L(/.)] , 



(2) 



(3) 



(4) 



appearing within the series coefficient of x^{iJi). The superscript Uf indicates the 
number of active flavours contributing to the evolution of the QCD coupling con- 
stant as and the running 6-quark mass. Thus, the series S[x, L] within (2) is of 
the form 



may be expressed in terms of the QCD expansion parameter 

x(/x)^a("/)(/x)/7r 
and n*'^-order polynomials of logarithms 

L{^^)^\og(\^^/J■''\^^)\'" 



S[x, L] = 1 + X (oo + oiL) + (bo + biL + h-iL^) 
+ x^ (co + CiL + caL^ + c^L^) + O (x^) . 



(5) 



The one- and two- loop coefficients oq, ai, 60, ^1, and 62 have been calculated for 
arbitrary n/ in [1]. For example, if n/ = 5, 



oo = 4.25360, ai = 5, 60 = 26.7848, W = 36.9902, 62 = 17.2917. 



(6) 



Note that the parameter within (3) and (4) is the (nonphysical) renormalization 
scale parameter. The full decay rate (2) must ultimately be independent of the 
choice of /x: 







dV \ 

P{x) 



[l-2^{x))-^ + f3ix)^ + 5^{x) 



(Pox^ + Pix' + P2X^ + ...), 



S[x,L], 



(7) 
(8) 



-f{x) = - (jox + 7ia^^ + 72^;^ + •••), (9) 

where, for Uf = 5, po = 23/12, /3i = 29/12, /?2 = 9769/3456, 70 = 1, 71 = 253/72, 
and 72 = 7.4195 [2]. If we expand the final hne of (7) perturbatively in powers of 
x''L"^, we find that 

= x{ai- 570) + x"^ {bi - aoPo - 5ao7o + 2ai7o - 571) + x'^L (262 - ao/?o - 5ai7o) 

+x^ (ci - 26o/3o - doPi - 56o7o + 26i7o - 5ao7i + 2ai7i - 572) 
+x^L (2c2 - 26i/?o - «i A - 5ai7i - 56i7o + 46270) 

+ x^L^ (3c3 - 262/3o - 56270) + 0{x^). (10) 

For Uf = 5 we see from (6) that the coefficients of x, x'^, and x^L vanish, as 
required. The requirement that 0{x^) coefficients in (10) also vanish may be used 
to determine the logarithmic coefficients ci, C2, and C3 within S[x,L\. For Uf — 5, 
we find from (6) and (10) that 

ci = 249.592, C2 = 178.755, C3 = 50.9144. (11) 

The coefficient cq, however, is RG-inaccessible to the order we are working, and 
must be obtained via a direct three-loop calculation. 

PADE-APPROXIMANT ESTIMATION OF 
NEXT-ORDER TERMS 

In the absence of a direct calculation, Pade-approximant methods provide a 
means for estimating next-order terms in a perturbative series, provided such 
a series obeys appropriate criteria generally believed to be applicable to field- 
theoretical perturbative series [4]. As an example of the power of such methods, 
consider the following perturbative series: 

This series is, in fact, the first five terms of the Maclaurin expansion of the function 
S{x) = sec(a;) — tan(a;). A [A^|M] Fade approximant to a series is the ratio of 
degree- iV and degree-M polynomials that replicates the known terms of a series. 
For example, the [2|2] approximant which replicates the series (12) is specified (up 
to a common scale factor for the numerator and denominator) by 



Specifically, the Maclaurin expansion of S'l^'^l is 

^[^1^1 = 1 - X + - + - + . . . . (14) 

The final series term in (14) may be regarded to be a Pade-approximant prediction 
for the next term in the series (12). This prediction is surprisingly accurate, as the 
underlying function's Maclaurin series expansion is 

115 2 

secfa;) — tanfa;) = 1 — x -\ — x^ H -x'^ x^ + . . . . (15) 

^ ^ ^ ^ 2 3 24 15 ^ ^ 

We see that the final series terms hsted in (14) and (15) are quite close: 19/144(= 
0.13194 . . .) ^ 2/15(= 0.13333 . . .). 

In the calculation of the B — > Xui~i^£ decay rate, the particular series with 
which we are concerned is of the form 

S[x, L] = l + Ri[L]x + R2[L]x^ + i^sW^;^ + ■ ■ ■ , (16) 

with 

Ri[L] ^ ao + aiL, R2[L] ^ bo + hL + b2L^ (17) 

fully determined by the previously-calculated coefficients (6), and with the next 
order coefficient 

i?3W = Co + CiL + 02^2 + 03^3 (18) 

assumed to be unknown. Of course, Ci, C2, and C3 are in fact accessible by RG 
methods [they are listed in (11)]; we will ultimately compare Pade-estimates to 
RG-determinations of these coefficients to ascertain the accuracy of the Pade- 
estimation procedure we utilise to obtain the RG-inaccessible constant cq. 

Wc seek a procedure to estimate R^ in terms of the known lower order co- 
efficients Ri and R2 appearing in (16). To do this, we assume that the er- 
ror of an [A^ — 1|1] approximant in estimating the N + 1*^ term in a series is 
inversely proportional to N. For example, the [0|1] approximant to the series 
1 + Rix + R2X^ + Rsx^ + ... is just 

1 — Rix 

The predicted coefficient of x^ in (19) is i?2°'^' = -^i) as opposed to the true 
coefficient R2. We find the relative error of the [0|1] approximant in predicting R2 
to be 

d[0|1] d f?2 d 

r _ -K2 - -Kl - -K2 

<^[o|i] = ^ = 5 ■ (20) 

-K2 -K2 

Similarly the [1|1] approximant to the same series is 

^ l + iR^ -R2/R.)x ^ 1 ^ ^ ^^^2 ^ /^2/^ A ^3 ^ _ _ _ _ (21) 

i — [K2/-tii)X ^ ' 



The predicted coefficient of is Rl/Ri, which we denote as -R3 and the relative 
error associated with this prediction is 

_ -R3 - -^3 RI/R1-R3 V 
0[i\i] = = . (22) 

If the relative error of an [A?^ — 1|1] approximant in predicting Rn+i is inversely 
proportional to N, i.e., if 

'^[^-^i^i = — r;;;-, — ^ "TV' ^^^^ 

as suggested in [5], we see from (20) and (23) that A = 1 — RI/R2. Since (23) also 
implies that = —A/2, we can use (22) to obtain an error-improved estimate 
of R3 entirely in terms of the known series terms Ri and R2 [6] : 

Rs = 2Rl/[Rl + RiR2]. (24) 



ESTIMATION OF THREE-LOOP COEFFICIENTS 

The problem with the prediction (24) for the calculation considered here is 
the incompatibility of R^ as a degree-3 polynomial in L (18) with the rational 
(non-polynomial) function of L one obtains via substitution of (17) into (24). One 
approach to this problem is to generate a least squares fit of (24) to the form 
(18) over the perturbative portion of the domain of !/(//), i.e., over the ultraviolet 
domain > mfe(/x). We define w = m^(/i)//i^, in which case L — — log(u') from 
(4), and obtain a least squares fit of (18) to (24) over the range < w < 1 [i.e., 
A* ^ ''^^(a*)] by optimizing the function [7] 

X^(co,Ci,C2,C3) 

[Cq - Ci log(w) + C2 log [w] - C3 log (w) j 



1 

dw 







R\{w) ^ R^{w)R2{w) 

(25) 

with respect to q, the polynomial coefficients of powers of L within (18). Upon 

substitution of the known [via (6) and (17)] functions R\{yS) = — ailog(w), 
i?2(w) =60 — 61 log(u)) + h^log^iyS) into the integrand of (25), we find from explicit 
numerical integration that 

X^(co, ci, C2, ca) = + 2c? + 24c^ + 720c^ + 2coCi + 4coC2 + 12coC3 



+12ciC2 + 48ciC3 + 24OC2C3 - 2237.27co - 5980.96ci 



- 24423.1C2 - 129077C3 + 6386542. (26) 

The four conditions dy^jdci = {i = {0,1,2,3}) provide the following Pade- 
estimates for q: 

Co = 198, ci = 261, C2 = 184, cg = 48.6 . (27) 

Note the excellent agreement of estimates for {ci, C2, C3} with the true values (11) 
obtained via RG methods. The relative errors Sci = (cj — cf^)/cf'-^ in the Pade 
estimates are respectively dci = —4.6%, 6c2 = 2.9%, Sc^ = 2.1%. These small 
errors suggest similar accuracy in the determination of the unknown coefficient Cq. 
Utilizing this value of cq in the L = limit, appropriate for the = mi,{ii) choice 
of renormalization scale in (1), we find that the series within (1) is now improved 
to 

t{b^ XJ-Ve) = K{4.17GeVf[l + 0.304 + 0.137 + 0.075] = K{l913GeV^). 

(28) 

The estimate (28) is obtained using the central value of Chetyrkin and Stein- 
hauser's [8] determination of mh{mi,) = 4.17±0.05 GeV an estimate which followed 
from a central value as{Mz) — 0.118 [also utilized to evolve x{ijl) in (28)]. 



RESIDUAL SCALE DEPENDENCE 

The scries within the rate (28) is evaluated at ^ = A. 11 GeV and appears to 
be the first four terms of a positive-term perturbative series. This positive-term 
character is evident for ji > A.YI GeV as well, the domain of L(//) for which L 
is positive, since the coefficients a^, 6^, and q (including the estimated value for 
Co) are all positive. If convergent, a truncated positive-term series necessarily 
underestimates the series sum. For values of /x larger than A.17GeV, this un- 
derestimation is even more severe, because of the positive growth of L{fj,). Thus 
r(// ^9 GeV) = K{1733GeV^) < r(// = 4.17 GeF), as evident from (28). 

As noted in [1], the perturbative series does not become locally flat until /i 
is well into the infrared region. This region [fi < 7725(7715)] can be accessed via 
threshold matching conditions [9] to a running coupling constant and running 
mass characterized by four active flavours: 



x^^^ [mb{mb)] 

X^^^ [7776(7776)] 1 + 0.1528 (x^^^ [m6(7776)])^ + 0.633 (x^^^ [7776 ( 7776)])' 

777^^^ [7775(7775)] 



777^^ [7775(7775)] 



1 + 0.2060 (x^^^ [7775(7775)])^ + 1.95 (x^^^ [7775(7775)])' 



(29) 



(30) 



If as{M,) = 0.118 and mf^(mfe) = 4.170 GeF, the conditions (29) and (30) imply 
that (4.170 Gey) = 0.07155, x(^) (4.170 GeF) = 0.07162, andmJ^^(4.170Gey) = 
4.177 Gey. Moreover, with 4 active flavours the two-loop constants characterising 
the rate (5) are no longer as given in (6); instead one finds for n/ = 4 that [1] 

bo = 25.7547, h = 38.3935, 62 = 17.7083, (31) 

and, from pcrturbative RG-invariance (10) with constants /?j and 7j characterized 
by four active flavours, 

ci = 263.84, C2 = 194.23, Cg = 54.109. (32) 

One can also obtain estimates of the Uf = 4 three- loop coefficients via optimization 
of the function (25): 

Co = 181.5, ci = 277.3, C2 = 197.6, C3 = 51.86. (33) 

The startlingly close agreement between (32) and (33) for the RG- accessible co- 
efficients ci, C2, C3 is indicative of comparable accuracy for the estimate of cq in 
(33). 

In rcf. [10], it is shown that the = 4 values for q are consistent with 
continuity of the rate T{B — > X^i'u^) across the Uf = 5 threshold at 4.17 GeV. 
The rate is shown to increase with decreasing /i until a local maximum is reached 
at = 1.775 Gey [10]: 

)="• 

r [// = 1775 GeV] 



= K 



m[^^ (1.775 Gey)l [1 - 0.6455 + 0.2477 - 0.0143] = ir(2071 GeV^). (35) 



The optimisation condition (34) ensures that = 1.775 GeV is the choice of 11 
for which there is minimal residual scale sensitivity of the rate, an established 
criterion for the choice of scale [3]. Note that the perturbative series within (35) is 
now an alternating series whose 3-loop term is much smaller than that in (28). The 
alternation of sign becomes possible because L{/j) is negative when fi < mh{fi), as 
evident from (4). If the series within (35) continues to be both alternating and 
montonically decreasing, the rate at /i — 1.775 GeV can be shown to be bounded 
from below and above: 

K X 2071 GeV^ < T [/x = 1.775 GeV] <K x 2122 GeV\ (36) 



Remarkably, the "minimal-sensitivity" point defined by (34) is quite close to the 
point at which the three-loop term vanishes entirely, a point at which three-loop 
estimation of the series, if asymptotic, would be most accurate. With RG- values 



(32) for ci, C2, C3 and the Pade- estimate cq — 188.5 for cq [corresponding to 
minimization of with respect to Cq after exphcit incorporation of RG- values of 
Ci, C2, C3 into (25)], the three loop term (18) is seen to vanish at /i = 1.835 GeV^, 
at which point r(/i = 1.835 GeV) = K x 2069 GeV^. The agreement between this 
value for the rate with the minimal-sensitivity estimate (35) is striking. 

\Vub\: THEORETICAL UNCERTAINTIES 

The minimal sensitivity rate (35) is subject to the following AT/K theoretical 
uncertainties: truncation error (t.e.) 

\AT/K\t.e. = [m^^^ (1.775)]^ (0.0143) = 51GeV^, (37) 

error in estimating cq (|Aco/co| ~ |Aci/cf*^| ~ 5%^, 

\AT/Kl, = 35GeV\ (38) 
uncertainty in q;,(M^)(= 0.119 ± 0.002 [2]), 

^^/KU = ( ^_!2S ) G^^' (39) 

(asymmetric bounds are because (35) is based upon ^^(Mj) = 0.118), uncertainty 
in mb{mb){^ 4.17 ±0.05 [8]), 

Ar/i^U, = ( +1^5 ) GeV, (40) 
and the asymmetric nonperturbative (n.p) contribution to the overall rate [11]: 

Ar/X|n.p. = ( Is? ) GeV. (41) 

Inclusive of all these contributions, the result (35) becomes 

r^Kx (2065 ± U%)GeV\ (42) 

Using (1) for the coefficient K, the branching ratio for inclusive charmless semilep- 
tonic decay is seen to be 

„ , \ , ^ iKbpGl [(2065± 14%)G'el/5] 

R^r{B^ XJ-V,) /T(B anythmg) = 4|_£ \ , ^^.^^^ (43) 

in which case 

iKfcl = (0.0949 ± 0.0066)i?^/l (44) 

The central value of (44) is in agreement with Uraltsev's relation between \Vub\ and 
the branching ratio [12], and is indicative of an overall ±7% theoretical uncertainty 
in this relation. 
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